o 
<n: 
^— > ■ 

O: 

x- 



d 

(N 




32nd International Cosmic Ray Conference, Beijing 201 1 



•^.e Rat tQ 

Multiwavelength Modelling of the Globular Cluster Terzan 5 

I. Busching 1 ' 2 , C. Venter 2 , A. Kopp 3 - 2 , O.C. de Jager 2 , & A.-C. Clapson 4 

1 Institut fur Theoretische Physik, Lehrstuhl IV: Weltraum- und Astrophysik, Ruhr-Universitdt Bochum, 44780 Bo chum, 
Germany 

^Centre for Space Research, North-West University, Potchefstroom Campus, Private Bag X6001, Potchefstroom, 2520, 
South Africa 

3 Institut fur Experimentelle und Angewandte Physik, Christian- Albrechts-Universitdt zu Kiel, Leibnizstrasse 11, 24118 
Kiel, Germany 

4 Max-Planck-Institut fur Kernphysik, PO Box 103980, 69029 Heidelberg, Germany 
ib@tp4.rub.de 

Abstract: Diffuse X-ray emission has recently been detected from the globular cluster (GC) Terzan 5 [8], extending 
out to ~ 2.5' from the cluster centre [6]. This emission may arise from synchrotron radiation (SR) by energetic leptons 
being injected into the cluster by the resident millisecond pulsar (MSP) population that interact with the cluster field. 
These leptons may also be reaccelerated in shocks created by collisions of pulsar winds [4], and may interact with bright 
starlight and cosmic microwave background (CMB) photons, yielding gamma rays at very high energies (VHE) through 
the inverse Compton (IC) process. In the GeV range, Fermi Large Area Telescope (LAT) has detected a population of 
GCs [2], very plausibly including Terzan 5, their spectral properties and energetics being consistent with cumulative 
magnetospheric emission from a population of MSPs. H.E.S.S. has furthermore detected a VHE excess in the direction 
of Terzan 5 [3]. One may derive constraints on the number of MSPs, iVtot, and the radial profiles of the GC B-field, 
stellar energy density, as well as the diffusion coefficient using the spatially-resolved X-ray, high-energy (HE), and VHE 
fluxes. If the Fermi LAT flux is due to magnetospheric processes, it will scale with the number of visible gamma- 
ray MSPs, iVvis- The HE spectrum therefore provides an independent way of constraining the number of MSPs (since 
-Vtot > A r v is). Consequently, the synthesis of available multiwavelength data presents a unique opportunity to constrain 
several parameters of the GC Terzan 5. 
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1 Introduction 

The first prediction [4] of high-energy (HE) and very-high- 
energy (VHE) fluxes from globular clusters (GC) invoked 
inverse Compton (IC) processes of relativistic leptons in- 
jected by the embedded population of around 100 millisec- 
ond pulsars (MSPs) interacting with the background radi- 
ation fields. The MSPs are sources of leptons with ener- 
gies up to a few TeV, which are accelerated in their mag- 
netospheres by very large electric fields (e.g., [5]). Addi- 
tionally, these leptons may be further accelerated in shocks 
resulting from colliding pulsar winds. Thus, [4] assumed 
mono-energetic as well as power-law injection spectra (cor- 
responding to the above two acceleration scenarios), and 
calculated the upscattering of stellar photons (coming from 
the many old-type stars in the cluster core) and cosmic 
microwave background (CMB) photons as these energetic 
particles diffuse outwards. Neglecting synchrotron radia- 
tion (SR) losses in their calculation, it was found that GCs 
should be detectable as unpulsed, point-like sources for 
Cherenkov telescopes, depending of the number of MSPs 



(iVtot) in the cluster, as well as the particle efficiency r] e 
(fraction of MSP spin-down power converted into parti- 
cles). 

The pulsed gamma-ray component coming from particles 
accelerated inside of the MSP magnetospheres was next 
calculated [17], while a similar calculation to that of [4] 
of the unpulsed IC component was performed [18] using 
an injection spectrum calculated from first principles and 
which is the result of acceleration and curvature radiation 
(CR) losses in the MSP magnetospheres, assuming no fur- 
ther particle acceleration. It was found that 47 Tucanae 
and Terzan 5 may be visible for H.E.S.S., depending on 
the assumed model parameters, particularly N to t and clus- 
ter B-field B (the model fixing r\ e to ~ 0.007 due to the 
magnetospheric acceleration process). 

Recently, H.E.S.S. published upper limits on the VHE 
gamma-ray emission from the GC 47 Tucanae, allowing 
us to infer N to t ~ 30 - 40 for B ~ 10^G, but N to t be- 
coming quite larger for B < 5 /iG or B > 30/iG. Next, 
Fermi LAT detected HE emission from 47 Tucanae, with 
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the spectrum being consistent with collective pulsed emis- 
sion from about 50 — 60 MSPs in the cluster ([18] in- 
ferred TVvis ~ 50). Since then, Fermi has detected sev- 
eral GCs [1, 2, 11, 14], very plausibly including Terzan 5. 
Following these detections, an IC scenario was considered 
to explain the HE fluxes seen by Fermi [7], as an alterna- 
tive to the usual CR assumption. For certain parameters, 
the Fermi flux may be reproduced, also predicting spec- 
tral components that should be visible in the VHE domain. 
However, such unpulsed IC components seem less domi- 
nant in the case of the cluster NGC 6624, where the HE 
emission seems to come almost exclusively from a single 
pulsar in the GC, PSR J 1823 -3021 A [13]. 

H.E.S.S. has just announced a VHE excess in the direction 
of, but offset from the centre of, Terzan 5 [3]. In addition, 
diffuse X-ray emission from Terzan 5 was measured [8], 
peaking at its centre but smoothly decreasing outwards, 
possibly resulting from SR. Also, several radio structures 
have been identified in the vicinity of Terzan 5, although 
no reliable estimate could be made for the radio index [6]. 
Lastly, new measurements of Terzan 5's distance [9, 15], 
core radius, half-mass radius, tidal radius, and total lumi- 
nosity [12] have become available. 

In view of all these observational developments and mod- 
elling efforts, the aim of this paper is to model the SR flux 
components expected from Terzan 5 using updated param- 
eters and refined model assumptions. A preliminary IC flux 
calculation has been presented elsewhere [20]. 

2 The Model 

Some key conclusions have been reached [4] regarding the 
modelling of GCs which may guide future efforts: 

• Pulsar winds mainly interact among themselves, as 
stellar winds are confined to only a small circumstel- 
lar region, as inferred from a simple colliding wind 
model. 

• Since the cluster B-field should be less than B ~ 
10~ 4 G, and given the typical energy densities of soft 
photon fields, leptons injected into the cluster mainly 
lose energy by the IC process, as the latter energy 
densities exceed the magnetic energy density. SR 
losses may therefore typically be neglected. 

• Acceleration of leptons are limited by an advection 
process along the surface of the shock for typical GC 
parameters (and not the shock structure, or IC and 
SR losses), yielding maximum lepton energies of ~ 
4-40 TeV. 

• IC radiation decreases rapidly when moving out 
from the GC centre. This should result in point-like 
gamma-ray sources for Cherenkov telescopes. 

In contrast to the last point above, an alternative IC cal- 
culation [7] solving a cosmic -ray diffusion equation (and 



using a slightly different stellar photon energy density pro- 
file) predicted that most of the HE radiation comes from a 
region beyond the GC core, implying that GCs should be 
extended sources. Also, we will include SR losses in our 
calculation below (for more details, see [16]). 

In this paper, we perform a similar calculation, but using 
updated structural parameters, a much larger bolometric lu- 
minosity, and a distance of d =5.9 kpc [9, 12, 15] to calcu- 
late the radiation and escape losses and resulting unpulsed 
fluxes, assuming Bohm diffusion and B-fields of B = 1 /xG 
and B = 10 yuG. We use an energy density profile similar 
to that given by [4], assuming bright starlight and CMB 
as soft photon targets with central energy densities of sev- 
eral hundred eV/cm 3 (for a temperature T = 4 500 K) 
and 0.27eV/cm 3 (for T = 2.76 K). Importantly, we also 
use a power-law injection spectrum, normalized to the total 
power output of the particles 

L e = N tot r] e (E Iot ), (1) 

assuming r] e = 0.01 and an average MSP spindown of 
{E I0 t) = 2 x 10 34 erg s _1 . This spectral shape will re- 
sult from reacceleration of particles in the cluster, after 
having escaped from the MSP magnetospheres (see Sec- 
tion 3). Since the details of this process are not clear, its 
actual shape and energy cutoffs will have to be constrained 
by data. 

3 Propagation of Leptons 

For our study, we assume the leptons accelerated by MSPs 
in the GC core to be reaccelerated at shocks in the CG core, 
leading to a power-law injection spectrum, and then propa- 
gate diffusively. As the MPS provide an almost steady par- 
ticle source, the distribution function /, depending on the 
spatial distance r from the GC core and energy E, satisfies 
the equation 

-S = V-(«(r,25)V/)-^(£7/). (2) 
Here we use the form of the spatial diffusion coefficient, 

K(r,E) = K (r){E/lTeV) a (3) 

that is inspired by the form of the Galactic diffusion co- 
efficient with a = 0.6 in plain diffusion models. E = 
E(r, E) is the rate of lepton energy losses due to IC and 
SR that depends strongly on the distance from the GC cen- 
tre due to the spatial variation of the radiation field; S is 
the source term. For this inital study, we keep k = k(E). 
As the level of turbulence inside the GC core will be higher 
than in the surrounding medium, we will allow for a ra- 
dial dependence of the diffusion coefficient n{r) in further 
studies. Eq. (2) is solved numerically. 

4 Results 

Firstly, we scaled the pulsed CR component [18, 20] to fit 
the Fermi LAT data [2]. This implies a number of visible 
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Figure 1 : Radial profile of the synchrotron flux at 1 keV for different diffusion coefficients inside the GC. Comparing our 
calculations to the diffuse X-ray profile measured by [8], we find a best-fit value of n = 2.5 kpc 2 Myr _1 at 1 TeV. 
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Figure 2: Sample sky maps of the synchrotron flux, 
centred on the GC centre for a diffusion coefficient of 
2.5 kpc 2 Myr _1 at 1 TeV. The left panel is for 1 1 cm, while 
the right panel is for 1 keV. 

MSPs 7V V ; S ps 60 ± 30. This is also consistent with the es- 
timate of 180±?g [2], and formally presents a lower limit 
for iVtot, since A?tot > N v - ls . However, the pair-starved 
model may overpredict the magnetospheric CR flux by a 
factor of a few, and furthermore may not be valid for all 
MSPs. Indeed, light curve modelling [19] has shown that 
only a small fraction of the current gamma-ray MSP popu- 
lation may be described using the pair-starved model. This 
argues that iV v i s may be even larger. For comparison, there 
are 34 radio-detected MSPs in this cluster 1 , while it is esti- 
mated that the total number should be ~ 60 - 200 [10]. 
Using the H.E.S.S. -measured spectrum in conjunction with 
the HE data, one may next constrain N to t and the cluster B- 



field. The Cherenkov Telescope Array (CTA) may provide 
more stringent future constraints on these parameters. 

We were able to constrain the diffusion coefficient k using 
the radial profile of the diffuse X-ray emission [8] . Figure 1 
indicates different radial profiles of the SR flux at 1 keV, 
calculated for different values of k, while the associated SR 
flux sky maps at 1 1 cm and 1 keV are shown in Figure 2 
for the best-fit value of n = 2.5 kpc 2 Myr _1 at 1 TeV. 

5 Extended, Offset VHE Excess? 

The scenario where MSPs are the main injectors of rela- 
tivistic particles into the GC implies that the resulting mul- 
tiwavelength flux due to SR and IC should be centred on 
the GC if most of the MSPs are located within the core 
radius. Although this scenario can explain the source ener- 
getics (for typical GC parameters), it is intriguing that the 
VHE excess seen by H.E.S.S. [3] appears to be extended 
and offset from the GC centre. There may be a few reasons 
for this: 

• There may be radio-faint and / or off-beam MSPs 
outside of the GC core that inject leptons, but have 
not been observed to produce pulsed emission. 

• A bulk transport of relativistic particles to another 
acceleration site may have occurred. 

• The VHE excess may be due to a see-through back- 
ground source. 



1 . http://www.naic.edu/~pfreire/GCpsr.html 
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• Some other sources inside the GC beside MSPs may 
be responsible for / contribute to the VHE flux. 

In the first case, it is expected that MSPs have very wide 
radio and gamma-ray beams (e.g., [19]), so they will only 
be invisible if they are intrinsically faint. Also, the cumu- 
lative pulsed emission consisting of single-MSP emission 
pulsed at different periods, as well as the typical short MSP 
periods, complicate blind searches for new ones. Second, 
a transport of particles is expected to leave a trail of ra- 
diation, although identification of such an emission struc- 
ture may be limited by the angular resolution of the tele- 
scope in the case of gamma rays. Third, it is quite improb- 
able that the VHE emission will coincide with Terzan 5 
by chance, but not impossible. Lastly, detailed modelling 
will be needed to demonstrate that the predicted spectral 
normalization and shape correctly reproduce the data in 
case other sources of HE leptons are invoked, as the MSP 
scenario has been reasonably successful at explaining the 
pulsed spectra of GCs. These tentative ideas need to be 
developed further to assess their viability. It should also 
be remembered that Terzan 5 may be considered a non- 
typical GC, showing signs of a complex stellar formation 
history [9]. 

6 Conclusions 

One may derive important constraints on the GC Terzan 5 
by using multiwavelength data. The combination of HE 
and VHE data will already constrain N tot and B, for a 
given energy density and diffusion coefficient, using fluxes 
averaged over the cluster's extent. However, a more refined 
model including particle transport and radiation as a func- 
tion of radius will allow the prediction of SR and ICS flux 
profiles, in addition to the total fluxes. 

Modelling the CR component [20], we found that iV v i s « 
60 ± 30, consistent with the estimates presented in [10, 2], 
and representing a lower limit for N tot . 

Assuming that the electrons injected by the embedded 
MSPs are reaccelerated at shocks in the GC core, and using 
a diffusion coefficient with the same energy dependence as 
in models of Galactic CR propagation, we find that a dif- 
fusion coefficient of 2.5kpc 2 Myr _1 at 1 TeV describes the 
radial profile of the X-ray flux as reported by [8] best (as- 
suming a cluster field of B = 10/iG). 

Future work will include attempts to constrain the radial 
profile of the diffusion coefficient, energy density, and clus- 
ter field as well as the number of host MSPs and the particle 
injection spectra using radio, X-ray, and gamma-ray data. 
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